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Abstract 

For the plane, sphere, and hyperbolic plane we consider the canonical invariant 
determinantal point processes Z p with intensity pdv, where v is the corresponding 
invariant measure. We show that as p — > oo, after centering, these processes converge to 
invariant H 1 noise. More precisely, for all functions / £ H l {v) r\L l {u) the distribution 
of J2z£Z p f( z ) ~ n I fd v converges to Gaussian with mean and variance 




1 Introduction 

Determinantal processes are point processes with a built-in pairwise repulsion. They were 
first considered by Macchi ^Uj as a model for fermions in quantum mechanics, and have 
since been understood to arise naturally in a number of contexts, from eigenvalues of random 
matrices to random spanning trees and non- intersecting paths, see [T[ |2J IH [7| |Hl El U21 UH1 UZj • 
A point process Z on C is determinantal if for disjoint sets Di, . . . , Dk we have 



E [TT#(AnZ)l = / det(K{ Xi , Xj )) dn(x 1 )---dn(x k ), (1) 

J\k V / l<i,j<k 
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for each k > 1. Here K(x,y) is a hermitian symmetric measurable function and p is some 
reference measure. The integrand is often called the joint intensity or correlation func- 
tion of the point process. 

Conversely, if such a function K defines a self-adjoint integral operator K, on L 2 (A,p) 
which is locally trace class with all eigenvalues lying in [0, 1], then there exists a point process 
satisfying (HJ). In this case we speak of the determinantal process (K,p). In a weak sense, 
the points repel one another because the determinant vanishes on the diagonal. 

The processes we consider are defined by two properties. First, they correspond to 
a projection K to a subspace of analytic functions with respect to a radially symmetric 
reference measure. Second, their distribution is invariant under the symmetries of their 
underlying space A. The latter is either the complex plane C, 2-sphere §, or hyperbolic 
plane U. We will think of them as a subsets of C, or strictly speaking C U {oo}, though for 
the sphere § usually there is no harm in ignoring the point at infinity. 

These properties are uniquely satisfied by a family of processes indexed by a single density 
parameter p > on each space, see Krishnapur ^TJ Theorem 3.0.5], who discovered several 
remarkable properties of these processes. 

Planar model. Here A = C and for any p > 0, consider the kernel 



K p (z, w) = e pzw with respect to dpJz) = —e p|z| dz. 



Here, as in the sequel, dz stands for Lebesgue area measure on C. 

Tk=0 

in L 2 (C, p p ), and so the above pair defines a determinantal process of infinitely many points 
in the complex plane, see, for example [Zj. 

Spherical model. The space A is § = C U oo, the two-sphere. Now p is integer valued, 
p = 1,2,... and the (K, p) pair reads 

K p (z,w) = [ k ) ( z ™) k = (! + zw) p -\ with respect to dp p (z) = ^(1 + \z\ 2 )~^dz. 

k=0 ^ ' 71 

Note that the reference measure p p is typically not a constant multiple of the invariant 
measure v. K p is a projection kernel, onto the orthonormal polynomials \J ( p ~ 1 ^z k , k = 
0, . . . , p — 1 in L 2 (S, p p ). In this case, p is really the total number of particles. 

Hyperbolic model. Take A = U, the unit disk, which we identify with the hyperbolic 
plane. Let 

K p (z,w) = -y+i w ^ n res P ec t to dp p {z) = — (1 — Iz^^dz, 



for any p > 0. As in the planar model, (K p , p p ) defines a determinantal process with infinitely 
many points; the ortho normal polynomials in L 2 (iJ, p p ) being J ( p ~j i k )z k for k — 1, 2, ... . 

The intensity measure in any determinantal processes is K(z, z)dp(z). In the above 
models we find that K p (z, z)dp p (z) = ^du(z), where v is the invariant measure on A, 
unique up to constant. Here we use 

dv c {z) = dz, du n (z) = n~m? - and du u( z ) = 7T~TT?v?^ ( 2 ) 

(1 + \z\ ) \V — \z\ ) 

The distribution of the above processes is invariant under symmetries of the respective A, 
i.e. linear fractional transformations of C preserving the measure v^. 

Our main theorem concerns the linear statistics for the point process. Recall that for 
/ : A — ► R, the ivariant measure v — v\ and the intrinsic gradient V t we have 

imiW) = / iv ji 3 ^, \\f\Ww = [ \f\du. 

J A J A 

We say / is in H 1 ^) or L 1 ^) if these corresponding norm is finite. 

Theorem 1. For either the planar, spherical or hyperbolic model, let f G H l (v) H L l (v). 
Then, as p — > oo ; the distribution of 



converges to a mean zero normal with variance ^\[f\(jj x t v y 

Note that for both the limiting variance and the shift to make sense it is necessary to 
have / G H 1 n L 1 , so the theorem holds for the most general test functions possible. 

The fact that the variance is of order 1 manifests the advertised repulsion. The ifx-norm 
is conformally invariant, so one may replace the intrinsic gradient and intrinsic measure by 
the planar gradient and Lebesgue measure for the embedding. 

This work is partially motivated by the recent results of Sodin and Tsilerson ^H] on 
the three canonical Gaussian analytic functions (GAFs) with zero sets invariant under the 
symmetries of the plane, sphere, and hyperbolic plane. These processes are also indexed by 
a density parameter, and jTHj establishes asymptotic normality for the corresponding linear 
statistics, with / e Cq. What is striking is that for GAFs the variance actually decays as 
the density tends to infinity: 



Var ^2 f( z ) ~ const x P L \\^\\V(uy 

zez p 



,-Hl A f 112 
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Thus, the zeros of typical GAFs are more orderly than their determinantal counterparts. 

The determinantal processes studied here, while attractive solely on the basis of their 
invariance, also arise as matrix models. The planar case is really just the "infinite dimensional 
Ginibre ensemble" . If A is an n x n matrix of iid standard complex Gaussians, then as n f oo 
the point process of A-eigenvalues converges to the planar model, and p here corresponds to 
scaling. As for the spherical model, Krishnapur rjl] has proved that it coincides with the 
eigenvalues of A~ X B, where A and B are independent p x p Ginibre matrices. Further, for 
integer p, Krishnapur provides strong evidence that the hyperbolic points have the same law 
as the singular points of A + zAy + z 2 A 2 H — • in \z\ < 1 with again A , At, . . . independent 
p x p Ginibre matrices. 

In all three cases, Krishnapur provides natural random analytic functions for which Z is 
the set of zeros. Using an integration by parts argument, Theorem^ can be interpreted to 
say that the log absolute value of these analytic functions converges to the Gaussian Free 
Field. See Section 3 in ^3] for this relation in the Ginibre ensemble. 

Theorem^ also identifies the present as a companion paper to ^1 which treats the lim- 
iting noise for the Ginibre eigenvalues. The eigenvalues also define a determinantal process 
in C, see jS]. In it is shown that, along with an if 1 -noise in the interior of U similar to 
above, there is an H l l 2 (d\}) noise component in the corresponding n j oo central limit theo- 
rem. The invariance and lack of boundary effects in the three models considered here makes 
for essentially different proofs that are shorter and rely less on combinatorial constructions. 
The only overlap is the cumulant formula which is the starting point for both proofs. 

The main Theorem ^ is proved in three steps. In Section |21 we establish some general 
conditions under which (smooth) linear statistics are asymptotically normal, without com- 
puting the asymptotic variance. For this, the fact that the kernels is an analytic projections 
and their specific decay properties are crucial. In Section El we check that these properties 
are satisfied by our models. In Section 0] we determine the asymptotic variance and extend 
the convergence to general test functions. 

2 General conditions for asymptotic normality 

Taking a broader perspective, this section shows that under certain conditions satisfied by 
the models we are considering, linear functionals are asymptotically normal. 

Let B be a compact subset of C. Consider the following set-up. K p : B 2 — > R is a set of 
kernels indexed by p, which ranges in an unbounded subset of the positive reals. The kernels 
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here are the Hermitian and are with respect to Lebesgue measure; more precisely, if 
K p denotes the kernels outlined above, and k = dfi p {z)/dz is the density of the reference 
measure, then here and below 



KJz,w) = K p (z,w) (k(z)k(j/)) 



1/2 



(3) 



Kernel Properties 

The eventual asymptotic normality rests on the following asymptotic properties of K p as 
p — > oo; all limit statements and o(-) notations refer to this limit. Throughout, c denotes a 
numerical constants which may change from line to line. 



Uniform bound (UB). It holds 



\K P \\ao '■= sup \K p (x,y)\ < cp. 

x,ydB 



L 1 bound (L 1 B). For x,y G B we have a bound 

\K p (x,y)\ < (p p (x-y) with H^l^ < c. 



(4) 



(5) 



Interaction decay (ID). The above bounding function satisfies 

IIMVji = 



(6) 



Limited local analytic projection (LLAP) property. Assume that B C C. Fix 
i?2 compact so that Bi C B°. For p — 0, 1, 2 we have 



Similarly, 



sup 



sup 



x p K p (x,z)- / K p (x,y)y p K p (y,z)dy 
Jb 



K p (x, z)z v 



K p {x,y) y p K p (y,z) dy 



B 



0(1) 



oil). 



(7) 



• Covariance (CO). For any function F with bounded third derivatives and compact 
support in the interior of B we have Cov p (d ZtZ F, zz) = o(l). 

Note of course that Gov p (f,g) indicates the covariance of J2 z ezf( z ) an d J2zez9( z ) ^ n 
the Xp-process. The main proposition of this section (proved as Proposition |HJ) is: 
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Proposition 2. Suppose that the above conditions are satisfied. Then, for the corresponding 
determinantal process any linear statistic f with compact support in the interior B° and 
bounded third derivatives is asymptotically normal. We also have convergence of all moments. 

That these conditions are satisfied by the planar, spherical and hyperbolic models is 
delayed to the next section. Here we provide a lemma that sheds more light as to how 
condition LLAP arises. 

Lemma 3 (Analytic projections restrict to LLAP). Let K p : S 2 — > K be a kernel for 
the projection to the space of all analytic functions on the open set S C C with respect to 
measures \L p . For compact B C S, let K p denote the restriction of K p (x,y)(p p (x)p p (y)) 1 ^ 2 
to B x B. If K p satisfies UB and ID, then K p satisfies the LLAP. 

Proof. Note that since for each z, the function y i— > yK p (y, z) is analytic, it follows from the 
analytic projection property that 



K p (x, y)y p K p (y, z) dy = x p K p (x 1 z). 
Thus, for J7J) it suffices to show that 

o(l), (9) 



sup 



K p (x } y)y p K p (y,z) dy 

S\B 



where recall B 2 C B°. Setting s = y — z, there is a polynomial q of degree p so that 
\y p \ < <?(l s l) f° r an choices of z G B 2 ,y G B c . Also, q(\s\) < c\s\ 3 as soon as \s\ bounded 
away from zero. So, for z a positive distance from S\B, we have that 

\y p K p (y,z)\ < c\s\ 3 <p p (s), 

and by UB, ID, the absolute value of the left hand side of (jHJ) is bounded above by 

cp J \s\ 3 (p p (s)ds = c x p x o(p _1 ) = o(l). 

The proof of (JHJ is identical since K p is hermitian symmetric. □ 



Cumulants 

Recall that for any random variable X, the cumulants Cunifc(X), k = 1,2,..., are the 
coefficients in the expansion of the logarithmic generating function, 

logE[e^] = ^^Cum fc (A), 
k=i 
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and X is Gaussian if and only of Cum k (X) = for all k > 3. In any determinantal process 
(Kp,fjLp), the cumulants of the linear statistic f(z k ) have the explicit form, 

Cum fc) „(/) = — 5^ I ; . . '. I ; / [II f( x i) klK p( x h x i+i) dx i--- dx m, (10) 



m=l feiH hfcm=fe i=l 

fei km>l 



where x m+ i = x\ is understood, the integral ranges over m copies of the full space (here 
B), and again we are absorbing the reference measure /i p into the K p kernel. The structure 
behind formula (|1U|) has been employed in the past to establish asymptotic normality for 
determinantal processes with various assumptions on the regularity of /. See in particular 
the pioneering work of Costin-Lebowitz [3| and the later papers of Soshnikov, ^H] and |19j . 
While going through cumulants, the method here is quite different. 
We define the multiple integrals: for / a function of Xi, . . . Xj., 



Kf{f) = J 



f(x 1 , , x h ) K p (xj, x j+1 )dn(xj) 
i=i 



'ir. 



(the indices are mod k), and, as another shorthand, if the fa are all functions of one variable, 
we set 

Kp{h, ■ ■ ■ , fk) = K^iMx,) ■ ■ ■ f k (x k )). (12) 

Note that the cumulant (|T0)) is just a weighted sum of terms of the form K p (gi, . . . , g m ), ob- 
tained by partitioning {1,2, ... ,k} into m parts I\, . . . , I m of sizes k\, . . . k m and setting t/j = 
Yljeii fr H ence ; we more generally seek conditions for the vanishing of Cum p (/x, . . . , f k ), 
defined in the obvious way, for k > 3. 

The first step is a collection of estimates on the integrals f|12j) . The most fundamental 
of these are Lemmas 8 and 9 below. The former allows one to reduce the dimension in 
certain instances; the latter allows for the replacement of the test functions fa by their cubic 
appr oximat ions . 

Lemma 4. Assumptions L X B and UB imply that 

\K P (fi,...j k )\ <c P \\ Ahum 

1=2 

Proof. The integral is bounded above by 

/ \f(x{)\ ip p (x 2 - xi) ■ ■ ■ ip p (x k - Xk-i) dxi ■ ■ -dx k x ||/ 2 ||oo ■ • ■ Halloo x sup \K p (x,y)\. 

J B h x,yeB 



00 

II 1 ,, 
I oo- 
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Changing variables, yi = Xi — for i > 2 allows the remaining integral to be bounded 
above by 



fe-i 



\fi(xi)\(f P (y 2 ) ■ ■■<p P (yk) dx x dy 2 ---dy k = ||/i||i||<p P ||i 

Bx(2B) k ~ 1 

and the claim follows from assumptions (L 1 B, UB). □ 

Lemma 5. Assume L l B, UB, ID and that /, are bounded for all i. If any fi and fj are 
supported on disjoint compact sets, then K p (fi, . . . , f k ) — > 0. 

Proof. Now use the bound 

K p (h, ...,f k )< c||^||oo J (f P {x2 - xi) ■ ■ • V P ( x k ~ %k-i) dxi--- dx k 

where the integral on the right is over the product of the supports of the fi to ft- By 
adjusting c, we may insert \xi — Xj\ 3 (which is bounded below on the domain of integration) 
to produce 

Kp(h, ■ ■ ■ , fk) < cWKpWo,, <f p (y 2 ) ■ ■ ■<fp(y k )\yi+i + . . . + yj\ 3 dx^y 2 ■ ■ ■ dy k 

J Bx(2B) k - 1 



' Bx(2B) k ~ 1 

< cll-gjoo ^ / (p p (y 2 )---(p p (y k )\ye\ 3 dxidy 2 ---dy k 

e=i+1 JBx(2B)k-l 

< c x p x o(p _1 ) = o(l), 

after a change of variables (yi = x% — Xi-\, i > 1) in line one. □ 

Lemmas |U and El were developed for the following purpose. 

Lemma 6. Assume L l B, UB, ID, and LLAP. Let fx, ■ ■ ■ , f k be bounded in B, with, for 
some i, fi supported in a compact B\ C B° , and, for some j ^ i, fj(z) = z p for z e B, with 
p e {0,1,2}. Then 

Kp{fx, ...J k ) = Kp(fi, • • • , fj-i x fj, f j+1 ...,f k )+ o(l). 
Similarly, for fj(z) = z p we have 

K p (h, ...,/*) = Kp(h> ■■■■> /i-i> fi x /j+i ■■■Jk)+ o(l). 

Proof. We prove the first claim with fj = z p ; the proof of the second claim is identical. By 
the cyclic nature of K p , we may assume % = 1. Fix a compact set B 2 such that B\ d B2 C 
B 2 C B°. By the disjoint union decomposition 

(BxB)\ {B 2 x B 2 ) = {{B\B 2 ) x B)U(B 2 x (B\ B 2 )) 



and Lemma 03 we have the following (restrictions are placed only on functions with indices 
adjacent to j): 

Kpifli ■ ■ ■ j fk) — Kp(fl, ■ ■ ■ ) /j-i1b 2 ) /j) /j+i1b 2) • • • ) fk 



< 



+ 



Kp{fl, • • • , fj-llB\B 2 ,fj, fj+1, ■ ■ ■ , fk) 

Kp(fl, • • ■ , fj-llB 2 ,fj, /j+i1b\B 2 ) ■ ■ - i fk) = 0(1). 



Also, 

Kpifli ■ ■ ■ 1 fj-l^-B 2 , /j; /j+llB 2 ) ■ ■ ■ : fk) ~ ^p(/l> • • • 5 fj-l^-B 2 x /j; fj+l^B 2 , ■ ■ ■ 1 fk) (13) 

i=l 



where 

S(x, 2) = Qf y)^ p (y, z) rfy - x p K p (x, 2)) l(x,ze B 2 ). 
But, by the LLAP assumption, we have that 

sup \S(x, z)\ = o(l), 

x,z£B 

and, after the familiar change of variables yi = X{ — for i ^ j + 1, the argument used 
in Lemma 0] yields that the difference (JT3j) converges to zero. An identical application of 
Lemma gives 



Kpifi, • • • > fj-i 1 B 2 x /?, /j+i1b 2) ■ ■ ■ ,fk) - K p {fx, . . . , x /j-, . . . , / fc ) = o(l), 
which concludes the proof. 



□ 



We close this subsection by showing that the assumed conditions enable one to Taylor 
expand inside the K p integrals. 

Lemma 7. Assume that UB, L X B and ID hold. Let fj, 1 < j < k have bounded third 
derivatives. Then we have 



K p (h, •••>/*) = ^ E k , 



ok 
P 



h(xi)^[ft\^i-^r mi 

i=2 



+ o(l), 



m=0 m2+-+"it=n 

where we use the standard tensor notation for the full first and second derivatives. 
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^ } (: Cl )(: Cj -: Cl )^) 



Proof. Starting at £ = 2, we will step-by-step replace fi(xx)f 2 (x 2 ) • • • fe( x i), ^ > 2 by an 
approximation g t of degree 2 at X\. 

2 r i 

gt(x 1 ;x 2 ,...,x t ) = ^2 S ® 

m=0 m2+...+m£=m |_ i=2 

Note that all the are bounded on For the step-by-step replacement procedure we need 
to bound de = ge-ife — gt- Towards this end, let 

ft = fe( x u^t) = fe( x i) + fe( x i)( x t - x i) + f"( x i)( x e ~ x i)® 2 - 

Certainly 

\ff(xi,X£) - fe(xe)\ < c\xi -x^, 
for C independent of £,Xi,Xg. Also, 

\di\ = \ge-ift - 9e\ < \ge-i\\ft ~ ft \ + \gt-ift ~ gel- 
Let yt — X£ — x\. Since ge~i is bounded, we have 

\9i-i\\fi-fi\<c(\v2\' + ... + \vi\*) (14) 

for the range of j/j. As gi is produced from gt-\J\ by dropping all terms that are of degree 
3 or 4 in the y^ there is a constant c such that 

\9t+i-9ift+i\ <c(M 3 + -.. + M 3 ) (15) 
on the range of the yi. (Any monomial in y^ of degree 3 or 4 and coefficient 1 is bounded 
above on the compact range by the right hand side for c large enough). 
Now we write the difference 

\K P [gi-iixi, x 2 , ■ ■ ■ , x i-i)fe{xi) • ■ ■ fk( x k)] 

- K p [g e {xi, x 2 ,..., xi)ft+i{xt+i) • ■ ■ fk( x k)] \ 
= K p [d e (xi, x 2 , . . . ,xi)ft+i(xt +1 ) ■ ■ -fk( x k)] ■ (16) 
By UB and Z^B (JTfij) is bounded above by 

c||ifp||oo / (pp(x 2 -x 1 )---(fp(x k -x k ^i)\d i (x 1 ; x 2 ,...,x £+1 )\dx 1 ---dx k 

< c\\ k p\\oo / p p (y 2 ) ■ ■ ■(p p (yk)(y2\yi\ 3 )dx 1 dy 2 - ■ -dy k , 

JBx(2B) k - 1 7 

with (|14|) and (j!5|) used in the second line. Again by L X B and ID, this in turn is upper 
bounded by 

clBlpgioo || v?p x \y\ 3 || a < cp x o(p" 1 ) = o(l) 
as required. □ 
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Proof of the proposition 

The above bounds on K p made use of UB, L X B , ID, and LLAP. If we add CO to the mix, 
the result is the following. 

Proposition 8. Assume that K p satisfies conditions UB, L X B, and ID. For k > 3, let 
fi, . . . fk, be of compact support and have continuous third derivatives. If in addition CO 
holds, then Cum p (/i, . . . , fa) — > 0. 

To prove this, remember that Cum p (/ 1 , . . . , /*.) is a weighted sum of terms in the form 
K p (gi, . . . , g m ) each gj being a product of the underlying /'s (here m <k). Lemma [7| gives 

+ o(l), (17) 

in which Sj = Zi — z±, and we use complex coordinates Si, s«. For example, 

9?\zi)s? 2 = %%9i{zM + d-Ag^i)^ + 2d z d s g i (z 1 )s i s i . 

Further, the /c-fold integrals on the right hand side of (J 17)) are all of the form K° h (h(zi)o~iO-j) 
where h(zi) is a C 3 compactly supported (in B) function of Z\ and a rj = s v ,s v or 1. Since 
functions of at most three of the Zi-s are present in this integrand, Lemma |H1 with p = 
allows us to reduce it to an at most threefold integral: 

kfih^z^a,) = Kf+VMzxfatTj) + o(l), 

where d = 0, 1 or 2 is the number of distinct variables in a^aj. 
Next, two applications of Lemma El gives 

kf(h{z x )s 2 ) = k p {h,z)-k p (hz,i) 

= k p (h,z)-k p (hz) + o(i) = (i). 

Similarly, all except the SiSj terms vanish. Even among those, only half of the terms with 
% 7^ j survive, depending on the order of conjugation. Again by successive applications of 
Lemma |H1 

kf{h{ Zl )s 2 s z ) = k; 2 (s 1 h(z 1 )s 2 ) + (i) = o(i), 

while 

k; 3 (h( Zl )s 2 s 3 ) = kf{h{ Zl )s 2 s 2 ) + o(i) = k p (h, zz) - k p (hzz) + 0(1). (is) 

n 



K P {g 



E E 

£=0 £ 2 +-+£ k =m 



T^om 



Apart from the error, the latter equals — Cov p (h, zz). 

Therefore, aside from a constant term, the only 0(1) contributions to the cumulant sum 
are of the form ([180. The possible choices of h are: with F — /1/2 • • • fk, 

d = F x {d Z)Z gi)/gi 

= F x (d z g i )(d z g j )/(g i g j ), 

and our full K p formula reduces to K p (-, zz) — K p (- x zz) applied to 

m 

i=2 2<i<j<rn 

Reverting back to the original test functions /1, / 2) • • • 1 fk, this is a weighted sum of the 
functions 

F u = Fx (d z - z f u )jf u 
F uv = Fx(d z f u )(d- z f v )/(f u f v ). 

Since Cum p (/i, . . . , fk) is symmetric under permutations of indices of the f^s, it suffices to 
show that the total weight over the cumulant sum for each one of the two types of terms F u 
and F uv vanishes. Also, as each g$ is a product of ki of the fi, then Gi is a sum of ki terms of 
type F u and ki(ki — 1) terms of type F uv . Similarly, when i 7^ j, Gij is a sum of kikj terms 
of type F uv . Thus, the total number of terms of each type is given by 

type F u : k 2 + . . . + k m 

m 

type 

F U v '■ ^ ki{ki 1) + ^ kikj, 

1=2 2<i<j<m 

while each type appears in the cumulant sum with a different coefficient. 

Finally we invoke property CO: K p (d Z)S F, zz) — K p ((d z>z -F) x zz) — > 0. Since 

d z ,,F = F x (J2( d ^fu)/fu + ^dMzWMifM"), 

after subtracting d z ^ z F/k from each term of type F u it can be replaced by k — 1 terms of 
type F uv with the opposite sign. Thus, our final count is 

typeF™ : J^( k i ~ kk i) + Yl kik r ( 19 ) 

i=2 2<i<j<m 

That is to say, each m < k term in the cumulant sum is the same constant multiple of (J19[) . 
That this vanishes for k > 3 when summed over the full cumulant expansion is the content 
of the next lemma. 
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Lemma 9. For each m > 1 let ip(k, m,ki, . . . , k m ) be a real-valued function. With 



it holds that 



T fc (*0 = E 



(-1) 



m— 1 



<p(fc,m, k 1 ,...,k m ) 



m * — ' fci!---fc m ! 

m=l fc 1 + ...+fc m =fc 

fcl,...,fcm>l 



2<«<i<m 



0. 



i=2 



/or a// A; > 3. 

Proof. First realize, if we denote </? = 1 by 1, then 

oo 

log(e*) = £ 



m=l 

oo 



m 



1! 2! 3! 



+ 



X 



k=l 



which explains why the k > 3 cumulant of any constant is zero. Now set 

f = f( X , y) = e* + X ye* = 1 + ^ + + f!<l±M + . . 

so that the coefficient of the y term in the y-power series expansion of log / is 



and similarly, 



1 d 2 



= ^T fc (A;i + ... + k m ) 
y=o k=i 



x 



y=o k=i 

In order to obtain the pure quadratic sums we set 



oo 



, \ x , / x , 2 xn 1 , a;(l + J/) , a; 2 (l + 2 2 y) rr 3 (l + 3 2 y) 



2! 



3! 



The coefficient of the y term in this power series expansion reads 



= J2r k (k 2 1 + ... + k 2 jx k . 

y=o k=i 
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These series produce the types of terms we are after up to the fact that our cumulant 
expressions do not have the first coeffiicient k\. To omit this, the above may be modified as 
in 



81 



d Qog/)(e*-l) 



dy f-l 

00 

J2r k (k 2 + ... + k m )x k , 



y _u » m= i 



m—l 



y=0 



k=l 



and 



s 2 



d {logg)(e x - 1) 



dy g-1 



y=o k=i 



2 )x k 

ml ^ 1 



and finally 



1 d 2 (log/)(e--r 



2 d*y f-l 



y=° k=i 

Now we have an easy proof of the claim. Simply note that the right hand side of (J2Uj) equals 
the coefficient of the x k term in 

d x 2 

-x —Si + S 2 + Sn = —\ 



the latter being a straightforward computation. 



□ 



3 Properties satisfied by our models 

We now verify that the conditions for asymptotic normality of Section El hold for the three 
invariant models. First the simple bounds UB, L X B and ID are checked. Note that for 
the planar and hyperbolic models, property LLAP follows from Lemma EJ both families of 
kernels are projections onto the space of analytic functions on the plane or unit disk. For 
the spherical model, in which the kernel projects onto a set of finite degree, LLAP requires 
separate proof. 

The first three properties follow from a similar pointwise bound for each of the kernels 
in question. Recall the definition Q of K p (z,w). 

Planar model. We have 

7T 
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and so 

\ K ( Z ,w)\ = £ e P(^)-(k| 2 +l»| 2 )/2) = P £ -p\z-w\ 2 /2_ ^ 
P ' 71 IT 

That is, (p p (z) = ^e~ p ' z ' 2 / 2 , and it is immediate that conditions UB, L 1 B, and ID are satisfied 
for all z, w G C. 

Hyperbolic model. Denote t = (l — \z\ 2 )(l — \w\ 2 ) and s = (1—zw) 2 . Then \s\ = t+\z— w\ 2 , 
and we have 



K p (z,w) = ^(t/s) 



7TS 



If \z\ < R < 1, then \s\ G [(1 — i?) 2 ,4], and, assuming p > 2, we get 



\K p (z,w)\ = -?-\t/s\^ = -B- 

7T\S\ 7l\s\ 



1 - 



\z — w\ 



< 



1 - 



\z — w\ 



p-1 



< 



p 



(22) 



-0>-l)|*-to| 2 /8 



(1-RY 



--■ <p P (z-w). 



Thus, conditions UB, Z^B, and ID are satisfied for all \z\ < R, \w\ < 1. 



Spherical model. Now put t = (1 + H 2 )(l + \w\ 2 ), s = (1 + zw) 2 . Then t — \s\ + \z — w\ 2 , 
and the kernel reads 

K p (z,w) = -?-(8/t)W 2 . 

Tit 

If, both \z\ and \w\ < R, then t G [1, 6] with b = (1 + -R 2 ) 2 , and 



|-Kp(*»«>) 



7Tf' ' ' 7rf 



1 - 



Z — W\ 



(23) 



as soon as p > 2. Just as before, we have UB, L 1 B, and ID when |z|, \w\ < R. 
Rounding out the basic properties we have: 

Lemma 10. The spherical model restricted to a ball {z < \B\} has the LLAP property. 

Proof. Fix p > 0, and assume that p > 1 + p. We consider a truncated kernel, which is a 
projection to the space of polynomials of degree at most p — 1 — p with respect to the same 
measure as K p , that is p p . That is, we introduce 



,£±1 
7TT 2 



p-l-p 



fc=0 



(zw) k , 
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with again t = (1 — |-2| 2 )(1 — |w| 2 ). This truncated kernel is shown to have LLAP, and then 
the truncation is shown to make a negligible difference. 
First note that 

/ K p (x, y)y p K p (y, z)dy = x p K p (x, z), (24) 

since K p is a projection into the space of polynomials of degree at most p — 1 and for z fixed 
and y p K p (y, z) is Ppiy) 1 ^ 2 times such a polynomial. Also, since for z fixed, K p (y, z) — K p (y, z) 
is orthogonal to all powers at least p — 1 — p of y on any radially symmetric set, we have 



y\<B 



K p (x, y)y p (K p (y, z) - K p (y, z))dy = 0. 



(25) 



Next, using the bound ( p fe 1 ) < { p ' 1 k ~ p )p p we find 



,e+l 

TTl 2 



P 



K P (z,w)\<p p (l + \zw\) 



p-i-p 



And, when \z\ < b < B < \w\, we have (1 + \zw\) 2 /t < a 2 for some a < 1. Thus, 



pP +1 



1 + \zw\ 



1 p-l-p 



t 1 ' 2 



Similarly, for \z\ < b < \B\ < \w\, we have 



\K p (z,w)\<^ 



1 + \zw\ 

tV2 



p-1 



< P ~ a o-\ 



By Hermitian symmetry this implies that 



sup 

\x\,\z\<b 



K p (x,y)y p K p (y,z)dy 



\B\<\y\ 

Finally, an easy estimate shows 



< p p+2 a 2p - 2 - p 



B\ 





y 


V 


(1 + 


y 


2)p/2+2 



dy = o{\). (26) 



sup \x p K p (x, z) — x p K p (x, z)\ = o(l), 

\x\,\z\<b 



(27) 



and (|24*j) . (|2l)|) . (|2l)j) . (|2*7j) together imply the first part of the LLAP property. The second part 
follows from Hermitian symmetry. □ 



4 Asymptotic variance and general test functions 

Our goal is to prove asymptotic normality with explicit variances for any / G L 1 R H 1 . We 
do this by proving normality and determining the variance asymptotics for an || ■ ||#i-dense 
set of functions and then giving a uniform variance bound for all functions. 
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First note the general formula valid for all bounded f,g of compact support: 

Cov p (f,g) = J f(z)g(z)K p (z,z)dz - J Jf(z)g(w)K p (z,w)K p (w,z)dzdw 

f{ z ){9{z) - g(w))\K p (z,w)\ 2 dzdw, 

which, after symmetrization, reads 

CoVp(f,9) = \ J J(f(z)-f(w))(g(z)-g(w))\K p (z,w)\ 2 dzdw. (28) 

Lemma 11. The subset of smooth functions with compact support not containing oo is 
|| • \\ H i-dense in H 1 ^) fl L}(v). 

Note that this subset is not dense in /J 1 (U), only among H 1 fl L 1 functions: harmonic 
functions h in U are H l -orthogonal to any compactly supported /. This can be seen via an 
integration by parts, moving the gradient from / to produce a Ah = 0. 

Proof. Replacing / by (/ A b) V (—6) and letting b — > oo shows that bounded functions 
are dense. Then convolving a bounded / with a smooth probability density approaching <5 
shows that bounded C 3 functions are dense. 

First consider the planar or hyperbolic case, and use the invariant gradient V t , measure 
v and distance dist t . We may apply a sequence of smooth cutoff functions g r to / which are 
equal to 1 on the ball of radius r but are compactly supported and have |V t g r | < 1. Let 
h r — 1 — g r . We have 

V t (/-/<7r)=/lrVJ + /VA. 

and therefore 

\\f-fg r \\ 2 H i < 2 J \h r VJ\ 2 + \fVX\ 2 dv{z) 

< 2/ \VJ\z)\ 2 dv{z)+2\\f\\ O0 f \f{z)\dv{z) 

J&ist, (0,z)>r «/dist t (0, z)>r 

these converge to for bounded / G L 1 fl H 1 as r — > 00. 

For the sphere, we again consider smooth /, and note that adding a constant to / does 
not change its if ^norm (formally, the space H 1 consists of equivalence classes of functions 
which differ by a constant). Adding a constant does not change the fact that / e Li either, 
as the invariant measure v§ is finite. So we may assume that /(oo) = 0, and by smoothness 
and compactness f(z) < Cf dist t (oo, z). We now take g £ (z) = ((dist t (z, oo)/e — 1) V 0) A 1, 
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which is supported on points at least e away from oo. Also, |V t <7(z)| < 1/e and vanishes for 
2 more than 2e away from oo. As before, we have 

f f c 2 f dist, (oo, z) 2 
Wf-f9e\\ 2 m < 2/ \VJ(z)\ 2 du(z)+2 -i ^-U-du(z) 

•/distj,(oo,z)<2e J6ist l (oo,z)<2E £ 

Both terms converge to when e — > 0, as required. □ 

Lemma 12 (Asymptotic variance for a dense set). Let f and g be C l and of compact 
support in A, where A = C for the plane or the sphere or U for the hyperbolic plane. Then 

Cov p (f,g) -> -^(f,9) Hl . (29) 

Proof. It suffices to compute lim^oo Var p (/) for / e Qj as the covariance may be identified 
by substituting f + g for /. First, by Taylor's theorem with remainder there is a bounded 
non-negative function e(r) tending to as r j for which 

Var„(/)- / / (V/(^)-(w-z)) 2 |A^, W )| 2 d^; 

< / / — w|)|z — u>| 2 2 (,2 — w) dzdw = o(l). 

Now examine the remaining integrand 

I(z,w) := |V/(z)| 2 |^ - W | 2 cos 2 (^)|i^ p (^, W )| 2 , 

where 8(z,w) is the angle between f(z) and w — z, under the change of variables w 
z + p~ l l 2 w' . Pointwise, in each of the three models, we have 

1_ _u/ 

where if)(z) = 1 (plane), = 1 + \z\ 2 (sphere), = 1 — \z\ 2 (hyperbolic plane). This would result 
in the limiting formula for the variance: with 9' denoting the limiting angle between z and 



1 

, 2 \K p (z,z+ ^-)\ 2 -+ ^77^exp(-K| 2 M,)). 



w', 



Var„(/) - ^ J^\Vf(z)\ 2 J^cos 2 (6') ^e-^ 2 /^dw'dz (30) 

i p poo p2n -i p 

= ^/iWWr/ / cos 2 (6')\wfe-^ 2 d6'd\w'\dz = — \Vf(z)\ 2 dz. 

ll{ J B JO JO 47r J B 

On the other hand, z and w ranging in a bounded set and ||V/||l°o < oo, the right hand 
side of 

\I(z,z + —)\ < cK|e-l™'l 2 / c 

is integrable on B x C (this again uses (|2"Tj). (j22j), and (J23j) ). Therefore, dominated conver- 
gence validates (f3T)J) and completes the proof. □ 



Corollary 13. The property CO holds for out models: Cov p (dzzf, zz1b(z)) — > for f £ C 3 ; 
supp(/) C -B° and B compact. 

Proof. Consider a smooth g satisfying 1#/ < g < 1^, where 5' is a neighborhood of supp(/). 
Lemma (|12j) with implies 

4nGoVp(d gg f,gzz) -> (d zz f,gzz) Hl = -(f,d zz (gzz)) Hl = 0. 

For the last equality, note that V/ and V<9 22 -(gzz) have disjoint support. 
Now 

Cov(<9 22 -/, (l fl - <?)) = K p (d z - Z f, (1 B - </)) - K p (d z - Z f x (l fl - #)), 

where the second term vanishes, and the first one converges to by Lemma and the fact 
that the arguments have disjoint support. □ 

Lemma 14 (Variance bound). Let f £ H l (v) flL 1 (z/). There is a universal c > so that 
for all p > 1 we have 

Var p (/) <c [ \VJ(z)\ 2 du(z). (31) 

J A 

Proof. By considering the negative and positive parts of / separately, we may assume / > 0. 
In each of the three models by the invariant version of (|28jl we have, for / bounded and 
compact support 

Var p (/) = i J J\f{z) - f(w)\ 2 \K L (z,w)\ 2 du(z)du(w) } (32) 

where K L (z,w) = K(z,w)(r](z)r](w))~ 1 ^ 2 and i] = du(z)/dz is the density of the invariant 
measure (J2J. 

Repeating the identities in (}2"%|) for the invariant K\ and using that for p fixed K L is 
bounded, we get that (f3~2"J) extends to all / £ L 2 (u), in particular for bounded / £ L l (y). 
Now replace the nonnegative / £ L l (y) by f n = f A n. Let V(/) denote the right hand 
side of (|32p. Since |/ n (^) — /n(w)| is monotone increasing in n, the monotone convergence 
theorem gives V(f n ) — > V(/). We also have Var p (/ n ) — > Var p (/): the mean converges to a 
finite limit as / £ L l (v) and the second moment converges by the monotone convergence 
theorem. Thus (|32|) holds for all / > in L 1 ^), although a priori both sides may be infinite. 

For isometries T of A (i.e. linear fractional transformations preserving v) we have 

\K\T(z),T(w))\ = \K'(z,w)\. (33) 

A simple way to check ()33|) is to write \K L (z,w)\ directly as a function of the single variable 
|T 2 (w)|, where T z is the isometry taking z to 0; such an expression is clearly invariant. It is 
also possible to get from the invariance of the process and the covariance formula (|28|). 
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Let 7 2U) be a geodesic connecting uj and z that proceeds at speed s zw = dist t (z, w) given 
by the invariant distance between z and w. Then we have 

2 



\f{z)-f{w)f 



< si 



\VJ\ 2 {l zw {t))dt. 



Thus Var p (/) is bounded above by |V t /| integrated against the measure 



dm 




A J A JO 



<W) (Q\ Kt ( z > w )\ 2dt du ( z ) rfz/ M- 



This is defined in an invariant way, so it must be a some a(p) £ (0, oo] times the invariant 
measure. In fact, 'it is the invariant convolution on the symmetric space A (see, for example 
jH] for background) of the standard invariant measure with the radially symmetric measure 



dm = \ 




A JO 



sLK,(t){C)\K"{0,w)\ 2 dt dv{w), 



and therefore 



a(p) = tf(A) = - / s 2 0w \K^O,w)\ 2 du(w) 



dt 



Ra 



7i / rsQ r \K L (0,r)\ 2 r)(r)dr 



7T 



77(0) 



Ra 



rs 2 0r \K(0, r)\ 2 dr. (34) 



Here R\ = 00 or 1 is the radius of the planar model for A. A direct computation shows that 
for all three models (|34p is bounded by an absolute constant, verifying our claim. □ 

Proof of Theorem^ Corollary (fT3*|) shows that condition CO holds, and the other condi- 
tions have been checked in Section El For / £ C 3 of compact support Proposition gives 
asymptotic normality and Lemma fT2l gives the limiting variance, so we have 



£,(/):=£/(*)-- / \f\d" A/Yo,-^- 



(35) 



Lemma El allows us to extend the preliminary conclusion (|33|) to the advertised result. 
For any / £ H 1 ^) (and of appropriate support) there is a sequence of f £ £ Cq with 
11/ ~~ feWn 1 — > as e — > 0. Moreover, Lemma IT41 implies that the family {Z p (f)} is tight and 



also that 



< E 



Z p (f)-Z p (f £ ) <c / |V t (/-/ ( 
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The right hand side can be made small at will. Now, choosing a subsequence p' over which 
Z p (f) has a limit in distribution, we find the Fourier transform of that limit is as close as we 
like to that of a mean zero Gaussian with variance J A |V t /| 2 . (The full limit for Z p (f £ ) 
exists for any e > 0). Since this appraisal is the same for any subsequence p', we have pinned 
down the unique distributional limit of Z p (f). □ 
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